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ABSTRACT 



We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Thanks to its great sensitivity, 
Planck is able, for the first time, to detect SZ emission up to r ~ 3 x Rsoo- We test previously proposed spherically symmetric models for the pressure 
distribution in clusters against the azimuthally averaged data. In particular, we find that the Arnaud et al. (2010) "universal" pressure profile does 
not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r < Rsoo'- by r = 2xRsqo it underestimates the 
QQ ■ observed y profile by a factor of a 2. This may indicate that at these larger radii either i) the cluster SZ emission is contaminated by unresolved SZ 
' sources along the line of sight or ii) the pressure profile of Coma is higher at r > /J500 than the mean pressure profile predicted by the simulations 
£Sj , used to constrain the models. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west 
7-H ■ and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western 
• • ■ feature was previously noticed in the ROSAT PSPC mosaic by Markevitch (2000) as well as in the radio. Using Planck y profiles extracted from 
e ' corresponding sectors we find pressure jumps of 4.5^ 5 and 4.9+ 7 , in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure 
■ jump conditions, we deduce that the shock waves should propagate with Mach number M w = 1.95+ ,- 4 , 5 and M se = 2.03^^ in the West and 
%| ■ Southeast, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales with small intrinsic 
' scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields 
fall off much more slowly with radius than previously thought. 

Key words. Cosmology: observations - Galaxies: cluster: general - Galaxies: clusters: intracluster medium - Cosmic background radiation, 
X-rays: galaxies: clusters 



1. Introduction 

The Coma cluster is the most spectacular Sunyaev-Zeldovich 
(SZ) source in the Planck sky. It is a low-redshift, massive, and 
hot cluster and is sufficiently extended that Planck can resolve 
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it well spatially. Its intracluster medium (ICM) was observed in 
SZ for the first time with the 5.5 m OVRO telescope (Herbig 
et al. 1992, 1995). Later, it was also observed with MSAM1 
(Silverberg et al. 1997), MITO (De Petris et al. 2002), VSA 
(Lancaster et al. 2005) and WMAP (Komatsu et al. 201 1) which 
detected the cluster with a signal-to-noise ratio of S /N - 3.6. 
As reported in the all-sky early Sunyaev-Zeldovich cluster pa- 
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per, Planck detected the Coma cluster with a S /N > 22 (Planck 
Collaboration 201 Id). 

Coma has also been extensively observed in the X-rays from 
the ROSAT all-sky survey and pointed observations (Briel et al. 
1992; White et al. 1993), as well as via a huge mosaic by XMM- 
Newton (e.g. Neumann et al. 2001, 2003; Schueckeret al. 2004). 
The X-ray emission reveals many spatial features indicating in- 
falling sub-clusters such as NGC4839 (Dow & White 1995; 
Vikhlinin et al. 1997; Neumann et al. 2001, 2003) , turbulence 
(e.g. Schuecker et al. 2004; Churazov et al. 2012) and further 
signs of accretion and strong dynamical activity. 

Moreover, the Coma cluster hosts a remarkable giant ra- 
dio halo extending over 1 Mpc, which traces the non-thermal 
emission from relativistic electrons and magnetic fields (e.g. 
Giovannini et al. 1993; Brown & Rudnick 2011). The radio 
halo's spectrum and extent require an ongoing, distributed mech- 
anism for acceleration of the relativistic electrons, since their ra- 
diative lifetimes against synchrotron and inverse Compton losses 
are short, even compared to their diffusion time across the clus- 
ter, (e.g. Sarazin 1999; Brunetti et al. 2001). The radio halo also 
appears to exhibit a shock front in the West, also seen in the X- 
ray image, and is connected at larger scales with a huge radio 
relic in the Southwest (Ensslin et al. 1998; Brown & Rudnick 
2011). 

In this paper we present a detailed radial and sector anal- 
ysis of the Coma cluster as observed by Planck. These results 
are compared with X-ray and radio observations obtained with 
XMM-Newton and the Westerbork Synthesis Radio Telescope. 

We use H = VOkms-'Mpc" 1 , Q„, = 0.3 and Q A = 0.7, 
which imply a linear scale of 27.7 kpc arcmin -1 at the distance 
of the Coma cluster (z = 0.023). All the maps are in Equatorial 
J2000 coordinates. 

2. The Planck frequency maps 

Planck (Tauber et al. 2010; Planck Collaboration 201 la) is the 
third-generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 2011) covers the 30, 44, and 70GHz bands 
with amplifiers cooled to 20 K. The High Frequency Instrument 
(HFI; Lamarre et al. 2010; Planck HFI Core Team 201 la) covers 
the 100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to 0. 1 K. Polarisation is measured in all but the 
highest two bands (Leahy et al. 2010; Rosset et al. 2010). A 
combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics (Planck Collaboration 2011b). Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
(Planck HFI Core Team 2011b; Zacchei et al. 2011). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration (201 lh-z). 

This paper is based on the Planck nominal survey of 14 
months. The whole sky has been covered two times. We refer 

1 Planck (http://www.esa.int/planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries: 
France and Italy) with contributions from NASA (USA), and telescope 
reflectors provided in a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 
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Fig. 1: Upper panel: Radial profile of y in a set of circular an- 
nuli centered on Coma. The blue curve is the best fitting sim- 
ple model to the profile over the radial range from 85 arcmin to 
300 arcmin. The model consists of a power-law plus a constant 
y ff. The best fitting value of y ff is shown with the dashed hor- 
izontal line. Two vertical lines indicate the range of radii used 
for fitting. Lower panel: The probability of finding an observed 
value of y — ycoma m a given annulus. The probability was esti- 
mated by measuring y in a set of annuli with random centers in 
any part of the image outside 5 R%)§. For the annulus between 2.6 
and 3.1 /?5oo (or 122 arcmin to 147 arcmin) the probability of get- 
ting ycoma by chance is ~ 3 10~ 3 (see Fig. 1). For smaller radii, 
the probability is much lower, while at larger radii the proba- 
bility of getting y in excess of ycoma is ~10% or higher. The 
results suggest that Planck detects the signal from Coma in nar- 
row annuli AR/R = 0.2 at least up to R max ~ 3 Rsqo- This is a 
conservative and model independent estimate. 



to Planck HFI Core Team (201 lb) and Zacchei et al. (201 1) for 
the generic scheme of time ordered information (TOI) process- 
ing and map making, as well as for the technical characteris- 
tics of the maps used. We adopt a circular Gaussian beam pat- 
tern for each frequency as described in these papers. We use the 
full-sky maps in the nine Planck frequency bands provided in 
HEALPix (Gorski et al. 2005) nside=2048 resolution. An er- 
ror map is associated with each frequency band and is obtained 
from the difference of the first half and second half of the Planck 
rings for a given position of the satellite. The resulting maps are 
basically free from astrophysical emission. However, they are a 
good representation of the statistical instrumental noise and sys- 
tematic error. Uncertainties in flux measurements due to beam 
corrections, map calibrations and uncertainties in bandpasses 
are expected to be small, as discussed extensively in Planck 
Collaboration (201 ld,c,e). 
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3. Reconstruction and analysis of the y map 

The Comptonization parameter y maps used in this work have 
been obtained using the MILCA (Modified Internal Linear 
Combination Algorithm) method (Hurier et al. 2010) on the 
Planck frequency maps from 100 GHz to 857 GHz in a re- 
gion centered on the Coma cluster. MILCA is a component sep- 
aration approach aimed at extracting a chosen component (in 
our case the thermal Sunyaev Zeldovich (tSZ) signal) from a 
multi-channel set of input maps. It is based mainly on the well 
known ILC approach (see for example Eriksen et al. 2004), 
which searches for the linear combination of the input maps that 
minimizes the variance of the final reconstructed map imposing 
spectral constraints. For this work, we apply MILCA using two 
constraints, the first to preserve the y signal and the second to 
remove CMB contamination in the final y map. Furthermore, we 
correct for the bias induced by the instrumental noise, and we si- 
multaneously use the extra degrees of freedom (dof) to minimize 
residuals from other components (2 dof) and from the instrumen- 
tal noise (2 dof). These would otherwise increase the variance of 
the final reconstructed y map. The noise covariance matrix is es- 
timated from jack-knife maps. To improve the efficiency of the 
algorithm we perform our separation independently on several 
bins in the spatial-frequency plane. The final y map has an effec- 
tive point spread function (PSF) with a resolution of 10' FWHM. 
Finally, to characterize the noise properties, such as correlation 
and inhomogeneities, we use jack-knife and redundancy maps 
for each frequency and apply the same linear transformation as 
used to compute the MILCA y map. The MILCA procedure pro- 
vides us with a data map y(k, m) (where k and m are the abscissa 
and ordinate, respectively) together with random realisations of 
an additive noise model dy(k, m), which is Gaussian, correlated, 
and may present some non-stationary behavior across the field 
of view. These maps are used to derive radial profiles and to per- 
form the image analysis as described below. 

We verified that the reconstruction methods GMCA (Bobin 
et al. 2008) and NILC (Remazeilles et al. 201 1) give results that 
are consistent within the errors with the MILCA method (see 
Planck Collaboration V 2012). 



3. 1 . Analysis of radial profiles 

In this paper, we present various radial profiles y(r) of the 2D dis- 
tribution of the Comptonisation parameter y(k, m). These allow 
us to study the underlying pressure distribution of the intraclus- 
ter medium of Coma. The y parameter is proportional to the gas 
pressure P = n e kT integrated along the line of sight: 

y = f w (i) 

m e c l J 

where n e and T are the gas electron density and temperature, cr r 
is the Thomson cross-section, k the Boltzmann constant, m e the 
mass of the electron and c the speed of light. All the radial pro- 
files y(r) are extracted from the y(k, m) map after masking out 
bright radio sources. In this work we model the observed y(r) 
projected profiles using the forward approach described in detail 
by e.g. Bourdin & Mazzotta (2008). We assume that the three- 
dimensional pressure profiles can be adequately represented by 
some analytic functions that have the freedom to describe a wide 
range of the possible profiles. The 3D model is projected along 
the line of sight assuming spherical symmetry and convolved 
with the Planck PSF to produce a projected model function 



f(r). Finally we fit f(r) to the data using a chi-square minimi- 
sation of its distance from the radial profiles y(r) + dy(r) de- 
rived from the MILCA map including its additive noise model 
y(k, m) + dy(k, m). This procedure allows us to decorrelate the 
additive noise fluctuations by projecting y{r) and the model f(r) 
onto the principal component basis of all realisations of the noise 
model, i.e. by diagonalizing its covariance matrix. 

The confidence intervals on each parameter are calculated 
using the percentile method; i.e., we rank the fitted values and 
select the value corresponding to the chosen percentile. Suppose 
we have 1000 realizations for a specific parameter £ which is al- 
ready ranked from bottom to top, the percentile confidence inter- 
val at 68% corresponds to [£i6o">> (my] - Notice that in this work 
the confidence intervals are reported with respect to the best fit 
value obtained by fitting the model to the initial data set. 

Note that, by design, the forward approach tests the capabil- 
ity of a specific functional to globally reproduce the observed 
data. For this reason, the error estimates represent the uncertain- 
ties on the parameters of the fitting function rather than the lo- 
cal uncertainties of the deprojected quantity. This technique has 
been fully tested on hydrodynamic simulations (e.g. Nagai et al. 
2007 ; Meneghetti et al. 20 1 0). 

3.2. Zero level of the y map and the maximum detection 
radius 

As a result of the extraction algorithm, Planck y maps are unde- 
termined to an arbitrary additive constant y ff. This constant is 
a free parameter in all our y-map models. This constant can be 
determined from the Planck patch by measuring the y value at 
very large radii where we expect to have small o no contribution 
to the signal from the Cluster itself. In particular, in the case of 
the 13.6° x 13.6° MILCA-based patch of the image centered on 
Coma, this constant is negative, as illustrated in Fig. 1. The ra- 
dial profile of y was extracted from the y map in a set of circular 
annuli centered at (RA J20 oo, DEC J2 ooo)= (12:59:47, +27:55:53). 
The errors assigned to the points are crudely estimated by calcu- 
lating the variance of the y map blocked to a pixel size much 
larger than the size of the Planck PSF. The variance is then 
rescaled for each annulus, assuming that the correlation of the 
noise can be neglected on these spatial scales. For a model con- 
sisting of a power-law plus a constant (over the radial range from 
85 arcmin to 300 arcmin) we find y off = -6.3 10~ 7 ± 0.9 10~ 7 . 
We note that the precise value of y^ depends weakly on the 
particular model used, and on the range of radii involved in the 
fitting. 

To determine the maximum radius at which Planck detects 
a significant excess of y compared to the rest of the image, 
we adopted the following procedure. For every annulus around 
Coma with measured y = ycoma we have calculated the distribu- 
tion of y = yrandom measured in A = 300 annuli of a similar size, 
but with the centers randomly placed in any part of the image 
outside the 5 x /?5oo circle around Coma, where R500 is the radius 
at witch the cluster density contrast is A = 500. When calculat- 
ing yrandom the parts of the annuli within 5 x R50Q were excluded. 
The comparison of ycoma with the distribution of y ra ndom is used 
to conservatively estimate the probability of getting y > ycoma 
by chance in an annulus of a given size at a random position 
in the image (see Fig. 1, lower panel). For the annulus between 
2.6 and 3.1 R500 (122 arcmin to 147 arcmin) the probability of 
getting ycoma by chance is w 3 10~ 3 (a crude estimate, given 
N = 300 random positions). For smaller radii the probability 
is much lower, while at larger radii the probability of getting y 
in excess of ycoma is ~10% or higher. We conclude that Planck 



3 



Planck Collaboration: The physics of the Coma cluster 




Fig. 2: The Planck y map of the Coma cluster obtained combining the HFI channels from 100 GHz to 857 GHz. North is up and west is to 
the right. The map is corrected for the additive constant y ff- The final map bin corresponds to FWHM = 10'. The image is about 130x130 
arcminutes 2 . The contour levels are logarithmically spaced by yfl (every 4 lines, y increases by a factor 2). The outermost contour corresponds to 
y = 2 x cr„ oise = 4.6 x 10~ 6 . The green circle indicates Rsm- White and black crosses indicate the position of the brightest galaxies in Coma. The 
white sectors indicate two regions where the y map shows a local steepening of the radial gradient (see Section 7 and Fig. 6). 



detects the signal from Coma in narrow annuli AR/R = 0.2 at 
least up to /? max ~ 3 x R500. This is a conservative and model- 
independent estimate. In the rest of the paper we use parametric 
models which cover the full range of radii to fully exploit Planck 
data even beyond /? max . 

4. XMM-Newton data analysis 

The XMM-Newton results presented in this paper have been de- 
rived from analysis of the mosaic obtained combining 27 XMM- 



Newton pointings of the Coma cluster available in the archive. 
The XMM-Newton data have been prepared and analyzed using 
the procedure described in detail in Bourdin et al. (2011), and 
Bourdin & Mazzotta (2008). We estimated the Y x - M gas x T 
parameter of Coma iteratively using the Yx - Msoo scaling rela- 
tion calibrated from hydrostatic mass estimates in a nearby clus- 
ter sample observed with XMM-Newton (Arnaud et al. 2010); 
we find R 500 m (47 ± l)arcmin w (1.31 ± 0.03) Mpc and we 
use this value throughout the paper. To study the surface bright- 
ness and temperature radial profiles we use the forward approach 
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described in Bourdin et al. (201 1) taking care to project the tem- 
perature profile using the formula appropriate for spectroscopy; 
i.e., we use the spectroscopic-like temperature introduced by 
Mazzotta et al. 2004. 



5. The Coma y maps 

The main goal of this paper is to present the radial and sectoral 
properties of the SZ signal from the Coma cluster. Here we de- 
scribe some general properties of the image; the full image anal- 
ysis will be presented in a forthcoming paper that will make use 
of all the Planck data, including the extended surveys. 

Fig. 2 shows the Planck y map of the Coma cluster obtained 
by combining the HFI channels from 100 GHz to 857 GHz. The 
effective point spread function (PSF) of this map corresponds to 
FWHM = 10' and its noise level is cr noiseI0 = 2.3 x 10~ 6 . 

To highlight the spatial structure of the y map, in Fig. 2 we 
overlay the contour levels of the y signal. We notice that at this 
resolution, the y signal observed by Planck traces the pressure 
distribution of the ICM up to /?5on. As is already known from 
X-ray observations (e.g. Briel et al. 1992; White et al. 1993; 
Neumann et al. 2003), the Planck y map shows that the gas in 
Coma is elongated towards the west and extends in the south- 
west direction toward the NGC4839 subgroup. Fig. 2 shows that 
the SZ signal from this subgroup is clearly detected by Planck 
(see the white cross to the southwest). 

The Planck y map of the Coma cluster also shows at least 
two marked local steepenings of the radial gradient, to the west 
and to the southeast with respect to the cluster center. For conve- 
nience, in Fig. 2 we outline these regions with white sectors. A 
steepening of the y radial gradient suggests the presence of a dis- 
continuity in the cluster pressure profile which may be produced 
by a thermal shock, as we discuss in Section 7. In Fig. 3 we 
show the Planck y map of the Coma cluster obtained by adding 
the 70 GHz channel of LFI to the HFI channels and smooth- 
ing to a lower resolution. The PSF of this map corresponds to 
FWHM = 30' which lowers the noise level by approximately 
one order of magnitude with respect to the 10' resolution map: 
0"noise3O = 3.35 x 10~ 7 . As for Fig. 2 the outermost contour level 
indicates y — 2 X a" no ; se 3o = 6.7 x 10~ 7 . Due to the larger smooth- 
ing, this map shows less structure in the cluster center, but clearly 
highlights that Planck can easily trace the pressure profile of the 
ICM well beyond R200 ~ 2 x Rsm (see the outermost circle in 
Figure 3). 



6. Azimuthally averaged profile 

Before studying the azimuthally averaged SZ profile of the 
Coma cluster in detail, we first show a very simple performance 
test. In Fig. 4 we compare the SZ effect toward the Coma clus- 
ter, in units of the Rayleigh- Jeans equivalent temperature, mea- 
sured by Planck and by WMAP using the optimal V and W bands 
(from Fig. 14 of Komatsu et al. 201 1). This figure shows that, in 
addition to its greatly improved angular resolution, Planck fre- 
quency coverage results in errors on the profile which are * 20 
times smaller than those from WMAP. Thanks to this higher sen- 
sitivity Planck allows us to study, for the first time, the SZ signal 
of the Coma Cluster to its very outermost regions. We do this 
by extracting the y radial profile in concentric annuli centered on 
the cluster centroid (RA,Dec)=(12:59:47,+27:55:53). 




Fig. 3: The Planck y map of the Coma cluster obtained by combining 
the 70 GHz channel of LFI and the HFI channels from 100 GHz to 
857 GHz. The map has been smoothed to have a PSF with FWHM = 
30' . The image is about 266x266 arcmin . The outermost contour cor- 
responds to y = 2 x o" no i S e3o = 6.7 x 10~ 7 . The green circles indicate Rsoo 
and 2 x R 5K> » R 20 o- 

We fit the observed y profile using the pressure formula pro- 
posed by Arnaud et al. (2010): 



(c 5m x)y[l + (c 5 o xrf-^ a, 

where, x = (R/R^qo). This is done by fixing /?goo at the best-fit 
value obtained from the X-ray analysis (R500 = 1.31Mpc, see 
section 4) and using three different combinations of parameters 
which we itemize below: 

- a "universal" pressure model (which we will refer to as 
Model A) for which we leave only Pq as a free parameter 
and fix c 500 = 1.177, y = 0.3081, a = 1.0510, /? = 5.4905 
(Arnaud et al. 2010); 

- a pressure profile appropriate for clusters with disturbed X- 
ray morphology (Model B) for which we leave Pq as a free 
parameter and fix c 500 = 1.08 3, y = 0.3798, a = 1.406, 
P = 5.4905 (Arnaud et al. 2010); 

- a modified pressure profile (Model C) for which we let all 
the parameters vary (except /?5oo)- 

The best-fit parameters, together with their 68% confidence 
level errors, are reported in Table 1 and the resulting best-fit 
models are overlaid in the upper left, upper right and lower left 
panels of Fig. 5, for models A, B, and C, respectively. We find 
that Eq. 2 fits the observed y profile only if all the parameters 
(except /?5oo) are left free to vary (i.e., Model C). 

We also fit the observed radial y profile using a fitting for- 
mula (Model D) derived from the density and temperature func- 
tionals introduced by Vikhlinin et al. (2006): 

P = n e xkT (3) 
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Fig. 4: Comparison of the radial profile of the SZ effect towards the 
Coma cluster, in units of the Rayleigh-Jeans equivalent temperature 
measured by Planck (crosses) with the one obtained by WMAP (open 
squares) using the optimal V and W band (from Fig. 14 of Komatsu 
et al. 2011). The plotted Planck errors are the square root of the diago- 
nal elements of the covariance matrix. Notice that profiles have been ex- 
tracted from SZ maps with 10' and 30' angular resolution from Planck 
and WMAP, respectively. 

where, 

2, x = 2 (rlr c T a 1 

n 2 

+ (4) 

and 

T(r) = T , jr . (5) 

°\\ + (rlr t n clb 

The best fit parameters, together with their 68% confidence 
level errors, are reported in Table 2. The resulting model is over- 
laid in the lower-right panel of Fig. 5. The above temperature 
and density functions contain many more free parameters than 
Eq. (2). All these parameters have been specifically introduced 
to adequately fit all the observed surface brightness and temper- 
ature profiles of X-ray clusters of galaxies. This function thus 
is capable, in principle, of providing a better fit to any observed 
SZ profile. Despite this, we find that compared with Model C, 
Model D does not significantly improve the quality of the fit. 

7. Pressure jumps 

Fig. 2 shows that Coma exhibits a marked local steepening of 
the radial gradient of its y signal at about 0.5 degrees to the west 
and a similar distance to the southeast of the cluster centre. This 
suggests the presence of discontinuities in the underlying cluster 
pressure profiles. To test this hypothesis and to try to estimate 
the amplitude and the position of the pressure jumps we use the 
following simplified approach: i) we select two sectors ii) we ex- 
tract the y profiles using circular annuli and iii) assuming spheri- 
cal symmetry, we fit them to a 3D pressure model with a pressure 



jump. The sectors are selected in order to follow, as closely as 
possible, the curvature of the y signal around the possible pres- 
sure jumps. This procedure is a bit arbitrary: the pressure jumps 
are not perfectly spherical symmetric, thus, the sector selection 
depends also on what is initially thought to be the leading edge 
of the underlying pressure jump. Choosing a sector that does not 
properly sample the pressure jump goes in the direction of mix- 
ing the signal from the pre and post pressure jump regions. The 
result is a smoother profile which, when fitted with the 3D pres- 
sure model, returns a smaller amplitude for the pressure jump 
itself. Thus, in the worst scenario, the measured pressure jumps 
represent a lower estimate of the jumps at the leading edges. In 
this work we set the centers and orientations of the west and the 
southeast sectors to the values reported in the first three columns 
of Table 3 and indicated in Fig. 2. In section 9.3 below we show 
that, within the selected sectors, the SZ and the X-ray analyses 
give consistent results. This indicates that i) these SZ-selected 
sectors are representative of the features under study and ii) that 
the hypothesis of spherical symmetry is a good approximation, 
at least within the selected sectors. 

We fit the profiles using a 3D pressure model composed of 
two power-laws with index rji and 772 and a jump by a factor Dj 
at the radius rj. It is important to note that, even if irrelevant for 
the estimate of the jump amplitude, the measure of both the slope 
772 and the absolute normalization of the 3D pressure at a given 
radius depends on the slope and extension of the ICM along the 
line of sight. To take this into account we assume that outside 
the fitting region (i.e. at r > r s , with r s = 2Mpc) the slope of the 
pressure profile follows the asymptotic average pressure profile 
corresponding to model C (i.e. 773 = f3 = 3.66; see section 6 and 
Table 1). The 3D pressure profile is thus given by: 

( DAr/rjy 1 ' r<rj 
P = P X { (r/rj)-K rj<r< r„ (6) 

{ (r s /rj)-^(r/r s )-^ r > r,. 

We project the above 3D pressure model integrating along the 
line of sight for r < lOMpc. 

The best fit parameters together with their 68% errors are 
reported in Table 3. Note, that the error bars on o are smaller 
than the angular resolution of Planck. As explained in the 
Appendix A, this is not surprising and is simply due to projection 
effects. 

In the left and right panels of Fig. 7 we show with a gray 
shadow the corresponding 3D pressure jump models with their 
errors for the west and southeast sectors, respectively. For con- 
venience in Fig. 6 we overlay the data points with the best fit 
projected y models and we show the ratio between the observed 
and the best fit model of the projected y profile in units of the 
relative error. This figure clearly shows that the pressure jump 
model provides a good fit to the observed profiles for both the 
west and southeast sectors. 

As reported in Table 3 the pressure jumps corresponding to 
the observed profiles are Dj = 4.5^^ and Dj = for the 

west and southeast sectors, respectively. 

8. SZ-Radio comparison 

In Fig. 8 we overlay the y contour levels from Fig. 2 with 
the 352 MHz Westerbork Synthesis Radio Telescope diffuse to- 
tal intensity image of the Coma cluster from Fig. 3 of Brown 
& Rudnick (2011). Most of the emission from compact radio 
sources both in and behind the cluster has been automatically 
subtracted. This image clearly shows a correlation between the 
diffuse radio emission and the y signal. 
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Fig. 5: Comparison between the azimuthally averaged y profile of the Coma cluster and various models. From left to right, top to bottom, we 
show the best-fit y models corresponding to the Arnaud et al. "universal" profile (A), the "universal" profile for merger systems (B), the modified 
"universal" profile (C, see 1), and the Vikhlinin et al. fitting formula (D, see. 2). For each panel we show in the Upper subpanel points indicating 
the Coma y profile extracted in circular annuli centered at (RAj2ooo,DECj2ooo)= (12:59:47, +27:55:53). The plotted errors are the square root of the 
diagonal elements of the covariance matrix. Continuous and dotted histograms are the best fit projected y model after and before the convolution 
with the Planck PSF, respectively. In the Lower subpanel we show the ratio between the observed and the best fit model of the projected y profile 
in units of the relative error. 



To provide a more quantitative comparison of the observed 
correlation, we convolve the the diffuse radio emission to 10 ar- 
cmin resolution to match the Planck y map. We then extract the 
radio and y signals from the innermost r < 50 arcmin region of 
the cluster and plot the results in Fig. 9. This is the first quan- 
titative surface-brightness comparison of radio and SZ bright- 
nesses 2 . We fit the data in the log-log plane using the Bayesian 
linear regression algorithm proposed by Kelly (2007) which ac- 



2 see Malu & Subrahmanyan (2012) for a morphological comparison 
between radio and SZ brightnesses in the Bullet cluster 



counts for errors in both X and Y. We find a quasi-linear relation 
between the radio emission and the y signal: 

y _ i f) (0.86±0.02) F (0.92±0.04) 
10-5 ~ iU r R ' 

where Fr is the radio brightness in Jy/beam (10 arcmin beam 
FWHM). 

Furthermore, using the same algorithm, we find that the in- 
trinsic scatter between the two observables is only cr scatter = 
(9.6 + 0.2)%. The quasi-linear relation between the radio emis- 
sion and y signal, and its small scatter, are also clear from the 
good match of the radio and y profiles shown in Fig. 10 obtained 
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Table 1: Best Fit Parameters for the Arnaud et al. (2010) pressure model (Eq. 4 and Eq. 5) 



model 
(Mpc) 


(io- 2 ) 


c 500 


y 


a- 


P 


^500 

(Mpc) 


A ("Universal") 
B ("Universal" merger) 
C ("Universal" all free) 


9 9^+0.03 
^'^ -0.04 

98+ 02 

"• 7o -0.02 
t ci +0.05 


1.17 
1.083 

9 qc+0.12 
z — 3 -0.01 


0.308 

0.3798 

<0.01 


1.051 
1.406 

i 91+O.O6 

^ -0.02 


5.4905 
5.49 

3 66 +0 - 50 


1.31 
1.31 
1.31 



Table 2: Best Fit Parameters for Model D pressure model (Vikhlinin et al. 2006) (eq. 3) 



Density Temperature 



n (10~ 3 cnr 3 ) 


1 


r (KeV) 


a 9+0.2 


r c (Mpc) 


0.18^1 


r,(Mpc) 


1 15 +01 


r,(Mpc) 


0-78^ 1 


a 


-0 20 +0 01 

" ,z -"-0.02 


a 


12 +002 1 
"• lz -o.oi I 


b 


9 +20 




59 +0 02 1 

"' J -0.02 1 


c 


+0 04 
"•"-o.oo 


7 


3 1 






e 


3 02 +012 1 






«02(cm" 3 ) 


fl | 







" The fit returns «02 = thus r C 2 and /?2 are arbitrary. 
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Fig. 6: Comparison between the projected y radial profile and the best fit y shock model of the west (left) and southeast (right) pressure jumps. 
Upper panel: The points indicate the Coma y profile extracted from the respective sectors whose centers and position angles are reported in table 3. 
The plotted errors are the square root of the diagonal elements of the covariance matrix. The continuous histogram is the best fit projected y model 
reported in Table 3. The two vertical lines mark the ±lcr position range of the jump. Lower panel: Ratio between the observed and the best fit 
model of the projected y profile in units of the relative error. 



by simply rescaling the 10' FWHM convolved radio profile by 
10° 86 x 10 5 . An approximate linear relationship between the ra- 
dio halo and SZ total powers for a sample of clusters was also 
found by Basu (2012), for the case that the signals are calculated 
over the volume of the radio halos. 

9. Discussion 

So far In this paper we have presented the data analysis of the 
Coma cluster observed in its SZ effect by the Planck satellite. 
In Section 5 and Section 6 we showed that, thanks to its great 
sensitivity, Planck is capable of detecting significant SZ emis- 



sion above the zero level of the y map up to at least 4 Mpc which 
corresponds to R » 3 x R^io. This allows, for the first time, the 
study of the ICM pressure distribution in the outermost clus- 
ter regions. Furthermore, we performed a comparison with radio 
synchrotron signals. Here we discuss our results in more detail. 



9. 1 . Global pressure profile 

To study the 3D pressure distribution of the ICM up to r — 3 - 
4 x /?5oo, we fit the observed y profile using four analytic models 
summarized in Table 1 and Table 2 (see Section 6). 
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Fig. 7: 68% confidence level range of the 3D-pressure model for the west (left panel) and southeast (right panel) sectors in Fig. 6. Gray shadows 
are the profiles derived from the Planck data. Red shadows are the profiles derived from the XMM-Newton data. 



Table 3: Best-fit parameters of the pressure jump model of Eq. (6) 



Sector 


"RA 
(J2000) 


"Dec 
(J2000) 


"Position angle 
(deg:deg) 


^0 

(10- 4 cnr 3 keV) 


(Mpc) 


Dj 


'/I 


m 


West 
Southeast 


13:00:25.6 
12:59:48.9 


+27:54:44.00 
+28:00:14.39 


340:364 
195:240 


8.8+°;i 
3.6!^ 


1 15+ 003 

1- -0.00 

92 +0 02 


4 5 +2 - 5 

A Q+0-7 


o o +0 - 20 

-0.00 
9 1+0.09 
^ -0.02 


1.28+™ 2 
1 00 +006 

lM -0.30 



" The RA and Dec indicate the center of curvature of the sectors from which the profiles have been extracted. 
* We fixed r s = 2 Mpc and 773 = 3.66 (see text). 



From the ratio plot shown in Fig. 5 we immediately see that 
the "universal" pressure profile (Model A) is too steep both in the 
cluster center and in the outskirts. The fit to the data thus results 
in an overestimation and underestimation of the observed SZ sig- 
nal at smaller and larger radii, respectively. The overestimation 
of the observed profile at lower radii is consistent with WMAP 
(Komatsu et al. 2011). This is expected since merging systems, 
such as Coma, have a flatter central pressure profile than the 
"universal" model (Arnaud et al. 2010). For merging systems, 
Model B should provide a better fit as it has been specifically 
calibrated, at r < R500, to reproduce the average X-ray profiles 
of such systems (Arnaud et al. 2010). Fig. 5 shows that this latter 
model indeed reproduces the data well at r < Rsoo- Nevertheless, 
as for Model A, it still underestimates the observed y signal at 
larger radii. The observed profile clearly requires a shallower 
pressure profile in the cluster outskirts as evident in Models 
C and D. This is important as the external pressure slopes of 
both Model A and B are tuned to reproduce the mean slope pre- 
dicted by the hydrodynamic simulations of Borgani et al. (2004); 
Nagai et al. (2007); Piffaretti & Valdarnini (2008) (from now on, 
B04+N07+P08). The Planck observation shows that the pres- 
sure slope for Coma is flatter than this value. This is also illus- 
trated in Fig. 1 1 where we report the pressure slope as a function 
of the radius in our models: we find that while at R — 3 x R^m 
the mean predicted pressure slope is > 4.5 for Models A and 
B, the observed pressure slope of Coma is » 3.5 as seen in 
Model C and Model D. In Fig. 12 we compare the scaled pres- 
sure profile of Coma with the pressure profiles derived from the 
B04+N07+P08 numerical simulations. We find that the Coma 



pressure profile at 2 x R^qq is already 2 times higher than the 
average profile predicted by these hydrodynamics simulations 
although still within the overall profile distribution which has a 
quite large scatter. The Coma pressure profile thus lies on the 
upper envelope of the pressure profile distribution derived from 
the above simulations. Note that the recent numerical simula- 
tions of Battaglia et al. (2011) predict a flatter mean pressure 
profile which appears to be more consistent with the Coma pro- 
file and, in general, with the Planck SZ pressure profile obtained 
by stacking 62 nearby massive clusters (Planck Collaboration V 
2012). It is beyond the scope of this paper to discuss in detail 
the comparison between theoretical predictions. Here we just 
stress that, at such large radii, there is the possibility that the 
observed SZ signal could be significantly contaminated by SZ 
sources along the line of sight. This signal could be generated 
by i) unresolved and undetected clusters and ii) hot-warm gas 
filaments. Contamination would produce an apparent flattening 
of the pressure profile. We tested for possible contamination by 
unresolved clusters by re-extracting the y profile excluding cir- 
cular regions of r — 5' centered on all NED identified clusters 
of galaxies present in the Coma cluster region. We find that the 
new y profile is consistent within the errors with the previous 
one, which implies that this kind of contamination is negligible 
in the Coma region. Thus, if there is SZ contamination it is prob- 
ably related to the filamentary structures surrounding the cluster. 
We note that from the re-analysis of the ROSATAll-Sky Survey, 
Bonamente et al. (2009) and Bonamente et al. (2003) report the 
detection of extended soft X-ray emission in the Coma cluster re- 
gion up to 5 Mpc from the cluster center. They propose that this 
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Fig. 8: Westerbork Synthesis Radio Telescope 352 MHz total intensity 
image of the Coma Cluster from Figure 3 of Brown & Rudnick (201 1) 
overlaid by the y contour levels from Fig. 2. Most of the radio flux 
from compact sources has been subtracted; the resolution is 133 arcsec 
x 68 arcsec at -1.5 degrees (W of N). The green circle indicates Rsw 
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Fig. 9: Scatter plot between the radio map after smoothing to FWHM = 
10' and the y signal for the Coma cluster. To make the plot clearer, we 
show errors only for some points. 



emission is related to filaments that converge toward Coma and 
is generated either by nonthermal radiation caused by accretion 
shocks or by thermal emission from the filaments themselves. 
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Fig. 10: Comparison of the y (black) and diffuse radio (red) global ra- 
dio profiles in Coma. The radio profile has been convolved to 10 arcmin 
resolution to match the Planck FWHM and simply rescaled by the mul- 
tiplication factor derived from the linear regression shown in Fig. 9. 
The radio errors are dominated by uncertainties in the zero level due to 
a weak bowl resulting from the lack of short interferometer spacings. 
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Fig. 1 1 : Comparison of the pressure slopes of the best-fit models shown 
in Fig. 5. The red, green, blue and grey lines correspond to Models A, 
B, C, and D, respectively. 

9.2. X-ray and SZ pressure profile comparison 

We can compare the 3D pressure profile derived from the SZ 
observations to those obtained by multiplying the 3D electron 
density and the gas temperature profiles derived from the data 
analysis of the XMM-Newton mosaic of Coma 

In Fig. 13 we compare the 3D X-ray pressure profile with the 
3D SZ profile of our reference Model C. We point the reader's 
attention to the very large dynamical range shown in the figure: 
the radius extends up to r — 4 Mpc probing approximately four 
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Fig. 12: Scaled Coma pressure profile with relative errors (black 
shadow) overplotted on the scaled pressure profiles derived from nu- 
merical simulations from Fig. 6 of Arnaud et al. (2010). 



orders of magnitude in pressure. In contrast, due to a combina- 
tion of i) relatively high background levels and ii) the available 
mosaic observations, XMM-Newton can probe the ICM pressure 
profile of Coma only up to a: 1 Mpc i.e., four times smaller radii 
than Planck, probing only m one order of magnitude in pressure. 

Due to the good statistics of both Planck and XMM-Newton 
data, we see that the pressure profile derived from Planck ap- 
pears significantly lower than that of XMM-Newton, even if they 
differ by only 10-15%. This low, but significant, discrepancy 
may be related to the fact that we are applying spherical mod- 
els to a cluster that has a much more complex morphology, with 
a number of substructures. A detailed structure analysis explor- 
ing these apparent pressure profile discrepancies is beyond the 
scope of this paper and will be presented in a forthcoming pa- 
per. Here we just show a comparison of the 3D pressure pro- 
files obtained from Planck and XMM-Newton in four 90° sectors 
centered on the cluster and oriented towards the four cardinal 
points (see Fig. 14). This shows that the pressure discrepancy 
depends strongly on the sector considered. In particular, we find 
that while in the north sector the Planck and XMM-Newton pres- 
sure profiles agree within the errors, in the west sector we find 
the largest discrepancy, up to 25-30%. As known from X-ray 
observations (see e.g., Neumann et al. 2003) the north sector is 
the one that is more regular, while the west sector is the one in 
which the ICM is strongly elongated with the presence of major 
structures. 
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Fig. 13: Comparison between the Planck and XMM-Newton -derived 
deprojected total pressure profiles. Upper panel: Blue line and light blue 
shadow: deprojected pressure profile, with its 68% confidence level er- 
rors, obtained from the X-ray analysis of the XMM-Newton data (see 
text). Black line and gray shadow: best fit and 68% confidence level 
errors from the Model C pressure profile resulting from the fit shown 
in Fig. 5. Lower panel: Ratio between the XMM-Newton and Planck - 
derived pressure profiles. The black line and the gray shadow indicate 
the best fit and the 68% confidence level errors, respectively. 



E 

o 
Q_ 

1.2 
0.8 



* 0.010 

E 

o 

Q_ 

1.2 
0.8 




100 1000 100 1000 

Rodius (kpc) Rodius (kpc) 



Fig. 14: Same as Fig. 13 but from profiles extracted in 4 90° sectors. 
From left to right, top to bottom we report the west (-45°, 45°), north 
(45°, 135°), east (135°, 225°)and south (225°, 3155°) sectors, respec- 
tively. 



9.3. Shocks 

In Section 7 we show that Coma exhibits a localized steepen- 
ing of its y profile in at least two directions, to the west and to 
the southeast. These suggest the presence of discontinuities in 
the underlying 3D pressure profile of the cluster. Using two sec- 
tors designed to follow the curvature of the y signal around the 



pressure jumps we estimate their amplitudes. This represents the 
first attempt to identify and estimate the amplitude of possible 
pressure jumps in the cluster atmosphere directly from the SZ 
signal. Interestingly, we find that similar features are observed at 
the same locations in the X-ray and radio bands. 

In Fig. 15 we compare the X-ray and radio cluster properties 
from the west and southeast sectors selected from the SZ image. 
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Fig. 15: Comparison of the X-ray and radio properties in the west (left column) and southeast (right column) sectors. Upper panel: Surface 
brightness and temperature profiles of the XMM-Newton mosaic. The continuous histograms show the best fit projected density and temperature 
discontinuous models to the data. The 3D pressure model is overplotted in Fig. 7. Lower panel: Radial profiles of 352 MHz radio emission 
at 2 arcmin resolution in the west (left) and southeast (right) sectors, after subtraction of radio emission from compact sources (see Brown & 
Rudnick 201 1). The two vertical lines mark the position range of the jumps. 



The X-ray surface brightness and temperature profiles have been 
derived from the XMM-Newton mosaic while the radio profile is 
extracted from the Westerbork observations at 352 MHz at 2 ar- 
cmin resolution. To guide the reader's eye, we mark, for each 
profile in the figure, the position of the pressure jump as derived 
from the analysis of the y signal (See Table 3). For both sectors 
we find that the X-ray surface brightness and radio profiles show 
relatively sharp features at the same position as the steepening 



of the Planck y profiles. This is also the case for the tempera- 
ture profile of the west sector. For the southeast sector, however, 
this evidence is less clear. Because it is located in a much lower 
signal-to-noise region of the cluster, the error of the outermost 
temperature bin is too large to be able to put a stringent con- 
straint on a possible temperature jump. 

To check if the X-ray features are also consistent with the 
hypothesized presence of a discontinuity in the cluster pressure 
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Table 4: Main parameters of the fit of the temperature and density models to the XMM-Newton data (see Eq. 8 and Eq. 8). The 
symbols M„, Mj, M n j, and M$z represent the Mach numbers derived from the X-ray density, temperature, and pressure (n x T), 
and SZ pressure jumps, respectively (see text). 



Sector 


r x 
(Mpc) 




M„ 


D T 




D„ x D T 


M„ T 




West 
Southeast 


1 17^+0.0003 

L.l '->_ 00 Q3 

0.9778+0 0002 


2 00+ 003 
7 4-1+0.02 

-0.02 


1 7^+0.03 
L ' '-'-0.03 

2 10+ 001 


3 0+ - 7 
1 3+ LS 

1 -'-0.6 


2 6+ ' 4 

^•"-0.4 

1 3+ 13 

1 -13 


6.0!!;t 
3-1-1.1 


9 q+0.2 
-0.2 
1 6+° 3 


1 95+ 045 

1-7 -0.02 

2 03+ 014 



profile we simultaneously fit the observed X-ray surface bright- 
ness and temperature profiles using the following discontinuous 
3D density and temperature models: 



n — no X 



and 



rx 

r > r x , 



r - t x, J D T(r/r x ) ([ r < r x 
1 - l0X Ur/r x )^ r>r x , 



(8) 



(9) 



where r x is the position of the X-ray jump and D„ and Dj are 
amplitudes of the density and temperature discontinuities, re- 
spectively. The above models are projected along the line of sight 
for r < 10 Mpc using a temperature function appropriate for 
spectroscopic data (Mazzotta et al. 2004). Notice that due to the 
poor statistics of the temperature in the southeast sector, for this 
profile we fix ^1 = £1 =0. This choice does not affect the deter- 
mination of the jump position r x which is mainly driven by the 
surface brightness rather than by the temperature profile. 

The best-fit position, density, and temperature jumps to- 
gether with their 68% confidence level errors are reported in 
Table 4. To make a direct comparison with the pressure jump 
measured from the SZ signal, in the same table we add the am- 
plitude of the X-ray pressure jump derived by multiplying the 
X-ray density and temperature models (P x = n e T). 

The best-fit surface brightness and projected temperature 
models are shown as histograms in Fig. 15. The best fit 3D 
P x model and its 68% confidence level errors are overplotted 
in Fig. 7. 

From Table 4 we see that the X-ray data from the west sector 
are consistent with the presence of a discontinuity both in the 3D 
density and 3D temperature profiles. Both jumps are detected at 
> 5<x confidence level and the pressure jumps derived from X- 
ray and from SZ are consistent within the 68% confidence level 
errors (Table 3 and Table 4 ). This agreement near the disconti- 
nuity is also illustrated by Fig. 7 which, in addition, shows that 
the 3D pressure profiles for the west sector derived from the SZ 
and the X-ray data are consistent not only near the jump but also 
over a much wider radial range. 

These results indicate that the feature seen by Planck is pro- 
duced by a shock induced by supersonic motions in the cluster's 
hot gas atmosphere. Assuming Rankine-Hugoniot pressure jump 
conditions across the fronts the discontinuity in the density, tem- 
perature and pressure profiles are uniquely linked to the shock 
Mach number. 

Table 4 shows that, consistent with the amplitude of the pres- 
sure jumps, the Mach number obtained from the SZ and X-ray 
pressure profiles are consistent within the lcr confidence level 
errors. Furthermore, the Mach number derived from the X-ray 
density and temperature profiles agree within the 2<x confidence 
level errors. This agreement supports the hypothesis that the 
west feature observed by Planck is a shock front. 



For the southeast sector Table 4 shows that the X-ray surface 
brightness profile is consistent with the presence of a significant 
discontinuity in the 3D density profile. Due to the low statis- 
tics, the temperature model returns large errors and Dj is not 
constrained (see Table 4). Thus, although consistent, we cannot 
confirm the presence of a temperature jump. Despite this we find 
that, as for the west sector, the pressure jumps and the pressure 
profiles derived from X-ray and from SZ are consistent within 
the 68% confidence level errors (see Fig. 7, Table 3, and Table 4). 
Finally, Table 4 shows that the Mach numbers derived from the 
amplitudes of the different 3D models are all consistent within 
the 68% uncertainties level. We would like to stress that this is 
true not only for Mj and M n j which, being directly connected to 
Dj, have relatively large errors, but also for M„ and M$z which 
do not depends on Dj at all. As for the west sector, this agree- 
ment supports the initial hypothesis that the southwest feature 
observed by Planck is also a shock front. 

Notice that the good agreement between the 3D pressure 
models derived from the X-ray and SZ data, both in the west 
and southeast sectors, indicate that, within the selected sectors, 
spherical symmetry is a good approximation to the underlying 
pressure distribution. 

We conclude this section by pointing the reader's attention 
to the fact that, even though the radio and X-ray observations 
have a much better PSF than Planck, Fig. 15 shows that the re- 
spective jumps in these observations appear smooth on a scale of 
=s 200kpc * 7'. As explained in detail in Appendix A this is sim- 
ply a projection effect. Despite its relatively large PSF, Planck is 
able to measure pressure jumps in the atmosphere of the Coma 
cluster as well as the other techniques. 

9.4. Quasi-linear SZ-radio relation 

In Section 8 we show that for the Coma cluster the radio bright- 
ness and y emission scale approximately linearly with a small 
scatter between the radio emission and thermal pressure. Due 
to the near-linear correlation, where line-of-sight projection ef- 
fects cancel out, we work here with volume-averaged emis- 
sivities. We first express the monochromatic radio emissivity 
[erg s" 1 cm" 3 str" 1 Hz" 1 ] as: 



e,. ~ n CRe B 



Qcrc 



B 



1 +«- 



(10) 



where a is the spectral index,B is the magnetic field, Z?cmb ~ 
3(1 + z) |i/G is the equivalent magnetic field of the CMB, and 
«CRe and QcRe are the density and injection rate of cosmic -ray 
electrons (CRe) , respectively. In general, <3cRe can be a function 
of position and electron energy, and will depend on the model 
of cosmic-ray acceleration assumed. In secondary (hadronic) 
acceleration models Dennison (1980); Vestrand (1982), the 
relativistic electrons are produced in collisions of long-lived 
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cosmic ray protons with the thermal electrons, resulting in 
QcRe n t hermai ncR P , where ncRp is the density of cosmic- 
ray protons. Recent models in this category (Keshet & Loeb 
2010; Keshet et al. 2010) require that, in contrast to nthermai, «CRp 
should be constant over the cluster volume in order to match the 
cluster radio brightness profiles. Pfrommer et al. (2008) show 
that there is strong CRp injection even in the cluster peripheries, 
due to the stronger shock waves there. Strong radio CRp diffu- 
sion and streaming within the ICM could also lead to a com- 
pletely flat CRp profile (EnBlin et al. 2011). In the limit where 
B » Z?cmb and assuming a « 1 (e.g. Giovannini et al. 1993; 
Deiss et al. 1997), this would lead to e r oc n t h e nnai K y/T. This 
is consistent with our observations ', especially since n me rmai( r ) 
varies much more than T(r) in the Coma cluster (see e.g. Arnaud 
et al. 2001; Snowden et al. 2008). Jeltema & Profumo (201 1) de- 
rive a lower limit for the average field in Coma of 1.7 p G, from 
limits on the Fermi y-ray flux. The y-ray analysis thus leaves 
open the question of whether Coma could be in the strong-field 
limit. 

However, the rotation measure observations of Bonafede 
et al. (2010) provide characteristic values of 4-5 pG for the 
combined contributions of the central diffuse cluster field and 
contributions local to each radio source (e.g. Guidetti et al. 201 1 ; 
Rudnick & Blundell 2003). The majority of Coma's volume, 
which is outside of the cluster core, is thus in the weak-field 
limit, which leads to 6, oc y B 2 /T. To remain consistent with 
the linear correlation found here, the magnetic field would 
thus need to be nearly independent of thermal density. The 
non-ideal MHD simulations of Bonafede et al. (2011) show a 
typical scaling of B oc t&fL,, which would yield e ( - oc y 2 - 2 /T 2 - 2 . 
This could make the secondary model inconsistent with the 
observations in the weak-field limit. 

Primary (re-)acceleration models assume that relativistic 
electrons are accelerated directly from shocks and/or turbu- 
lence generated in cluster mergers. The turbulent re-acceleration 
model (Schlickeiser et al. 1987; Brunetti et al. 2001; Petrosian 
2001) leads to a scaling of e r oc n me rmai T L5 oc y Vr in the B < 
Bcmb limit (Cassano et al. 2007) if one assumes B 2 oc n t hermai 
which is close to the simulation scaling results of Bonafede 
et al. (2011). Such a scaling relation is consistent with the ob- 
served correlation. However, in order to connect the CRe den- 
sity to n t hermai, primary models depend on a large number of free 
parameters, which are generally fit to match the observations. 
Recent attempts to reduce the number of assumptions by intro- 
ducing secondary CRe* as seed particles (Brunetti & Lazarian 
201 1, , see above) fail to reproduce the linear correlation in the 
weak-field limit. This is another manifestation of the problem all 
simple models have in explaining the large extent of cluster ra- 
dio profiles when compared to the X-rays and inferred magnetic 
fields (e.g., Dolag & EnBlin 2000; Govoni et al. 2001; Donnert 
et al. 2010; Brown & Rudnick 201 1). In future, robust measure- 
ments of the cluster's magnetic field profile, coupled with high- 
resolution radio/X-ray/SZ correlations, will be needed to rule out 
these naive models. 



3 In the case that a = 1 + 6, the relationship would be e, oc (y/T) 1+6 ' 2 , 
if we assume B oc yn thcrraal . E.g., if we use a = 1.2 (Giovannini et al. 
1993), we would expect y cc Te°- 9i , which is approximately our mea- 
sured value. However, we continue to use the term "linear relationship", 
with the understanding that the difference between our measured value 
and linear is consistent for our simple assumption about the spectral 
index. 



9.5. Pressure jumps and radio emission 

Shocks play an important role in the production of radio emis- 
sion. We expect that shocks created during cluster mergers will 
compress magnetic fields and accelerate relativistic particles. 
However, the radiating electrons will quickly lose their energy 
post-shock, and may not be visible for more than ~100 kpc 
behind the shock (e.g., Markevitch et al. 2005) given char- 
acteristic shock velocities and magnetic fields at //G levels. 
These shock-accelerated electrons, in shock-compressed mag- 
netic fields, have been proposed as the explanation for the ob- 
served polarized radio synchrotron radiation from cluster pe- 
ripheral relic sources (Ensslin et al. 1998). Lower fields do not 
increase the electron lifetimes, and can even decrease them at 
fixed observing frequency, because of inverse Compton losses 
against the CMB. Recent simulations show that the presence of 
cluster-wide turbulence following a major merger is maintained 
for a few Gyr at a few percent thermal pressure (e.g., Dolag et al. 
2005; Vazza et al. 2006; Kang et al. 2007; Paul et al. 201 1). This 
turbulence can re-accelerate mildly relativistic seed electrons, 
and is potentially responsible for the large-scale halo emission 
(see above). In addition, an extensive population of low Mach 
number shocks is also seen in simulations (e.g., Miniati et al. 
2000; Pfrommer et al. 2006) and could play an important role in 
particle re-acceleration. 

Shocks will also induce turbulence in the post-shock re- 
gion (~200-300 kpc). There are hints from the small-scale X- 
ray residuals (Fig. 3 of Schuecker et al. (2004)) that such tur- 
bulence may exist interior to the possible shocks seen in the 
west and southeast. For the western region, the combination of 
SZ/X-ray pressure jump, X-ray suggested turbulence, and ex- 
cess synchrotron emission, points toward a connection between 
turbulence and diffuse synchrotron emission. The details of that 
connection, however, are not clear. In addition to direct acceler- 
ation by turbulence, the post-shock synchrotron emission could 
be a result of a population of weaker, as yet undetected shocks, 
or freshly accelerated CRp interacting with the ICM in a region 
where turbulence has amplified the magnetic field (e.g., Dolag 
et al. 2005; Ryu et al. 2008; Kushnir et al. 2009; Keshet 2010). 
Synchrotron spectral indices and magnetic field measurements, 
in combination with reliable measurements of weaker shocks 
and turbulence, would be needed to discriminate between po- 
tential models. 

The expected rapid loss of radio emissivity post-shock 
can also help us understand why shocks are sometimes easily 
detected in the radio, and sometimes not. In the clearest cases, 
radio shocks are seen beyond any central halo as relatively thin 
structures known as peripheral relics (van Weeren et al. 2010) 
where they can accelerate relativistic electrons. Radio shocks 
may also be found at or near the edge of the halo, and would 
be characterized by a sharp, but low contrast, rise in brightness, 
while the post-shock emission blends in with the halo instead of 
falling off. The western shock described here in Coma, as well 
as suggested shocks at the edges of halos in Abell clusters 521 
and 754 (Giacintucci et al. 2008; Markevitch 2010; Macario 
et al. 2011) are likely examples of this case. Contrast effects 
will camouflage the appearance of shocks that are projected 
against any radio halo emission. This is probably the case for 
the y shock in the southeast, where the radio halo extends far 
beyond the shock. The Coma cluster thus hosts all three types 
of "radio shocks": a) the classic peripheral relic at a distance 
of 1.7 Mpc from the center (which Ensslin et al. 1998 and 
Brown & Rudnick 2011 suggest is an "infall" shock), b) the 
western shock at the edge of the halo, and c) the southeast shock 
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projected against the fading radio halo. 



10. Conclusions 

We present the SZ observations of the Coma cluster based on the 
Planck nominal survey of 14 months. The excellent sensitivity of 
Planck allows, for the first time, the detection of SZ emission out 
to at least R w 3 x . We limit our investigation to the radial 
and sectoral properties of the intracluster medium, and we study 
the pressure distribution to the outermost cluster regions. Our 
three main results can be summarized as follows: 

- the Coma pressure profile is flatter than the mean pressure 
profile predicted by the B04+N07+P08 numerical simula- 
tions and lies on the upper envelope of this specific simu- 
lated profile distribution. This trend has also been found in 
the pressure profile derived by stacking 62 nearby clusters 
of galaxies observed with Planck (Planck Collaboration V 
2012). 

- Planck detects a localised steepening of the y profile about 
half a degree to the west and also to the southeast of the 
cluster centre. Features in the X-ray and radio synchrotron 
profiles at similar locations suggest the presence of shock 
waves that propagate with Mach number M w = 1.95^9 ™ and 
M se = 2.03*qH in the west and southeast directions, respec- 
tively. 

- the y and radio-synchrotron signals are quasi-linearly cor- 
related on Mpc-scales with only small intrinsic scatter. 
This implies either that, unlike the thermal plasma, the en- 
ergy density of cosmic -ray electrons is relatively constant 
throughout the cluster, or that the magnetic fields fall off 
much more slowly with radius than previously thought. We 
detect a correspondence between the western y shock and 
previously reported radio/X-ray edge, and we argue that ei- 
ther the magnetic fields are strong in the cluster outskirts, 
which would permit the hadronic model to explain the ra- 
dio emission, or some sort of re-acceleration by turbulence 
or additional shock waves must operate in the region behind 
the detected outer shock structures. 

Even though this analysis is based on only about half of the data 
collected by Planck, our results represent a substantial step for- 
ward in the study of the physics of the Coma cluster. The full 
set of data collected by Planck, will not only improve the signal- 
to-noise by another factor ~ V(2) but also significantly improve 
our understanding of instrumental effects. Thus, we will be able 
to generate more accurate y maps, and more thoroughly unveil 
Coma's two-dimensional SZ structure and its filamentary envi- 
ronment. 
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Appendix A: Pressure profile discontinuities as 
seen by Planck 

In this section we show that the 10 arcmin angular resolution of 
Planck is still sufficient to detect and measure 3D pressure jump 
features in the Coma cluster. This is because at the cluster 
redshift, 10 arcmin ss 280 kpc, which is of the order of the 
smoothing induced by projection. To show this, we assume a 
discontinuous 3D pressure profile described by Eq. 6. We 
project this profile along the line of sight and we calculate the 
corresponding y profile as observed with an instrument with 

- i) infinite angular resolution; 

- ii) a 10' FWHM angular resolution, as for Planck. 

The results of this exercise are illustrated in Fig. A.l where we 
compare four different cases with and without pressure jumps. 
In the upper section of each panel of Fig. A. 1 we show the input 
3D pressure profiles. In the middle sections we show the 
projected pressure profiles without any smoothing (as black 
histograms). The lower sections show deviations from a single 
power-law fit. The panel on the left shows that, due to simple 
projection effects, the y profile appears smoothed with an 
equivalent smoothing scale of w 200 - 250 kpc. For Coma, this 
corresponds to an angular scale of 7 - 9 arcmin. In the same 
panel we overplot, as a red histogram, the y profile convolved 
with the Planck PSF This illustrates that the effect of the PSF 
smoothing is secondary with respect to projection effects. This 
indicates that there is only a modest gain, from the detection 
point of view, in observing this specific feature using an 
instrument with a much better angular resolution than Planck. 
Notice that the fact that the projection of the 3D pressure 
distribution onto a plane converts the sharp jump into a curved 
surface brightness profile allows us to recover the position of rj 
with an accuracy higher than the angular resolution of Planck. 
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Fig. A .2: Comparison of the 3D pressure model (green line) and the 
corresponding projected y profile, smoothed with the 10' beam. The 
dotted black line shows the expected y profile due to the inner power 
component, while the dashed line shows the contribution of the outer 
power-law component. The black solid line is the sum of these two com- 
ponents. For comparison, the blue line shows the same model, but not 
convolved with the 10' beam. In this plot Tj ~ 30'. Due to projection 
effects, the range of radii affected by the value of rj is of order rj itself. 
Since rj exceeds 10', many independent data points with large signal- 
to-noise ratio contributes to the determination of rj, allowing rj to be 
estimated with an uncertainty below the nominal angular resolution of 
the telescope. 

Fig. A. 2 clearly shows that the range of radii affected by rj 
value is of the order of rj itself. Since rj > 10', many 
independent data points with large signal-to-noise ratio are 
contributing to the determination of rj, driving the uncertainty 
well below the nominal angular resolution of the telescope. Of 
course the value of rj is still subject to systematic uncertainties, 
e.g. from our assumption of spherical symmetry. 
We use the above exercise to illustrate two practical ways to 
identify the presence of a possible underlying 3D discontinuous 
pressure profile, hidden behind some observed projected profile 
extracted in a specific cluster sector. 

The first way is to search for the actual pressure jump in the 
observed pressure profile. As projection smooths the profile, 
this needs to be done by looking at the profile extremes. We 
first notice that the outermost bins of the profile are practically 
unaffected by PSF smoothing. This is clear in the middle 
sections of Fig. A. 1 where the red and back histograms are 
similar in the outermost 3-4 bins. At this point, if we fit a line to 
the observed 3-4 outermost bins, this will give an indication of 
the un-convolved y profile slope at large radii. If we extrapolate 
this line to the center, we can easily highlight the presence of a 
pressure jump. This procedure is shown in the middle section of 
Fig. A. 1 with black straight lines. From Fig. A. 1 we can see 



16 



Planck Collaboration: The physics of the Coma cluster 




800 1104 1522 2100 800 1104 1522 2100 800 1104 1522 2100 800 1104 1522 2100 

R (kpc) R (kpc) R (kpc) R (kpc) 



Fig. A . 1 : Illustration of the effect induced by projection and by the Planck PSF on the y radial profile produced by an underlying broken power- 
law pressure profile with and without a pressure discontinuity. In this figure we fix Pq = 10. x l(T 4 cirr 3 keV, rj = 1.1 Mpc and, from left to right 
we consider four different cases: i) r]\ = , ni = 2, Dj = 4; ii) 771 = , ni = 2, Dj = 1; iii) rj[ = 2 , rj2 = 0, Dj = 4; iv) r/i = 2 , 772 = 0, 
Dj = 1. Upper Panels: The underlying 3D pressure profile; Middle Panels: The black and red histograms are the projected y profiles observed by 
an instrument with infinite angular resolution and with a PSF of 10 arcmin, respectively. The red line represents the best fit of a simple power-law 
to the entire convolved profile (red histogram). The black line is the same as the red line but considering only the three outermost projected profile 
bins. Lower Panels: Ratio between the PSF convolved and projected y profile and its best-fit power-law model (red histogram and lines in middle 
panel) in units of a relative error which, for this illustration, we set to 10%. 
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Fig. A .3: Results of the fit of the y profile extracted from three cluster sectors with a simple power-law. Left Panel) sector with the same angular 
size and extension as the west shock but pointing to the north, where there are no visible shock features; Middle Panel) profile of the west shock; 
Right Panel ) profile of the southeast shock. In the lower panel of each figure we report the ratio between the observed and the best fit model of the 
projected y profile in units of the relative error. The figure clearly illustrates that while the power-law gives a good fit for the north sector where 
no shock is present, it returns a poor fit in the west and southeast sectors. These two cluster regions require a discontinuity in the pressure jump as 
shown in Fig. 6. 



that, this procedure highlights the intensity variation due to a 
pressure jump (see first and third panels). In the case where we 
have no pressure jump (second and forth panels) the best-fit line 
to the outermost bins tends to closely follow the entire profile. 

The second way to highlight the presence of a pressure jump is 
to fit a line to the entire observed profile and to examine the 
residuals. This procedure is illustrated in the middle and lower 
sections of Fig. A. 1 . The red lines in the middle section are the 
best fit power-law relations to the entire observed profile (red 
histogram). The crosses in the lower panels indicate the 
differences, in units of the relative errors, between the PSF 
convolved projected y profile and its best fit power-law model. 



This figure shows that a 3D pressure jump induces a 
characteristic signature in the residuals. 

In Fig. A. 3 we apply this second technique to the Coma cluster 
by showing the fit of the y profile extracted from three cluster 
sectors with a simple power-law. The figure clearly illustrates 
that while the power-law gives a good fit for the north sector 
where no shock is present, it returns a poor fit in the west and 
southeast sectors. These two cluster regions require a 
discontinuity in the pressure jump as shown in Fig. 6. 
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